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Stellar evolution with rotation and magnetic fields: 
III: The interplay of circulation and dynamo 
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Abstract. We examine the effects of the magnetic field created by the Tayler-Spruit dynamo in differentially 
rotating stars. Magnetic fields of the order of a few 10* G are present through most of the stellar envelope, with 
the exception of the outer layers. The diffusion coefficient for the transport of angular momentum is very large 
and it imposes nearly solid body rotation during the MS phase. In turn, solid body rotation drives meridional 
circulation currents which are much faster than usual and leads to much larger diffusion coefficients than the 
magnetic diffusivity for the chemical species. The consequence is that the interplay of the thermal and magnetic 
instabilities favours the chemical transport of elements, while there would be no transport in models with magnetic 
field only. We also discuss the effects on the stellar interior, lifetimes and HR diagram. 
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1. Introduction 

Stellar rotation influences all the outputs of stellar evolu- 
tion and nucleosynthesis and several grids of models have 
been made for massive stars (Langer et al. 119991 Heger et 
al. l^DDUl Meynet & Maeder 12000 ) . The effects of rotation 
are even more important at metallicities Z lower than so- 
lar (Maeder & Mevnet |2U(J1|I . However, we do not know 
the role of the magnetic field in stellar evolution and it is 
the purpose of this series of works to approach this mat- 
ter, which is considered to be critical one (cf. Roxburgh 
I2003f) . Spruit (2002J proposed a dynamo mechanism op- 
erating in stellar radiative layers in differential rotation. 
This dynamo is based on the conjectured Tayler instabil- 
ity, which is apparently the first one to occur in a radiative 
zone (Tavler lWTHI Pitts & Tavler lT^5S| . For now, there is 
no empirical or observational proof of the existence of this 
instability. According to Pitts & Tayler H1986fl and Spruit 
H2002|l . even a very weak horizontal magnetic field is sub- 
ject to this instability, which then creates a vertical field 
component, which is wound up by differential rotation. 
As a result, the field lines become progressively closer and 
denser and thus a strong horizontal field is created at the 
energy expense of differential rotation. 

In a first paper (Maeder & Mevnet 1^0031 paper I), we 
have shown that in a rotating star a magnetic field can 
be created during MS evolution by the Spruit dynamo. 
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We have examined the timescale for the field creation, its 
amplitude and the related diffusion coefficients. The clear 
result is that magnetic field and its effects are quite im- 
portant. In the second paper (Maeder & Meynet 120041 
paper II) , a generalisation of the equations of the dynamo 
was developed. The solutions fully agree with Spruit's so- 
lution in the two limiting cases considered f Spruit I!j002fl . 
i.e. "Case 0" when the /i-gradient dominates and "Case 1" 
when the T-gradient dominates with large non-adiabatic 
effects. Our more general solution encompasses all cases 
of /i~ and r~gradients, as well as all cases from the fully 
adiabatic to non-adiabatic solutions. 

Paper II suggested that there is a complex feedback be- 
tween the magnetic instability, which generates the field, 
and the thermal instability which drives the meridional 
circulation (Maeder & Mevnet r2004|l . However, it was be- 
yond the scope of paper II to make numerical models of 
the interaction between circulation currents and dynamo. 
In this paper we account for this feedback, the main steps 
of which are the following: 

— Differential rotation creates the magnetic field. 

— Magnetic field tends to suppress differential rotation. 

— The absence of differential rotation strongly enhances 
meridional circulation, 

— Meridional circulation increases differential rotation. 

— Differential rotation creates the magnetic field (the 
loop is closed). 
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Thus, a delicate balance between the thermal and mag- 
netic effects occurs during evolution. The balance of these 
effects influences the transport of chemical elements to the 
stellar surface and the internal transport of angular mo- 
mentum, which determines the evolution of the surface 
rotation. 

In Sect. 2, we collect in a consistent way the basic 
equations of the dynamo. In Sect. 3, we give the basic 
expressions for the transport coefficients. In Sect. 4, we 
calculate numerical models for the interaction of the dy- 
namo and circulation on the internal stellar structure. In 
Sect. 5, the effects at the stellar surface and in particular 
the evolution of abundances are analysed. Sect. 6 gives the 
conclusions. 

2. The general equations for the dynamo 



with a frequency equal to the Brunt-Vaisala frequency N. 
In a medium with both thermal and magnetic diffusivity 
K and rj, the Brunt-Vaisala frequency becomes (Maeder 
& Meynet mnH) 



v/K + 2 



N' + Nf, 



(2) 



with the radiative diffusivity K = .f'^'^K^ and with the 
following partial frequencies 



-^V^ and = ^(Vad-V) 



(3) 



The present set of equations for the dynamo based on 
Tayler-Spruit instability has three important advantages 
with respect to the system of equations by Spruit ( 2002 ). ~ f^pPN 



The thermodynamic coefficients S and (f are S = 

- (Ilf#)p,^ ^^^'^ = (fTf^)T,p- is the pressure scale 
height. The restoring oscillations will have an average den- 
sity of kinetic energy 



(4) 



1. The equations by Spruit apply to the non~adiabatic 
case. As most of the stellar interior is adiabatic, this 
situation is not appropriate. The present equations 
fully treat the problem whether it is adiabatic or non- 
adiabatic. 

2. Either the /i-gradient was considered to dominate or it 
was absent. Such a situation is not adequate when mild 
mixing effects, which produce small /i-gradients, are 
studied. In addition, as always in all problems of mix- 
ing, the /i-gradients need to be very carefully treated, 
otherwise the results are incorrect. 

3. There is no need to interpolate between the two asymp- 
totic solutions called "Case 0" or "Case 1". The in- 
terpolation was made by Spruit l|20021) as a function 
of the parameter q, which defines the amount of dif- 
ferential rotation (see definition given in Eq. ^ . This 
is hazardous, as it expresses the interpolation of ef- 
fects concerning the /i-gradient and the degree of non- 
adiabaticity as a function of differential rotation. This 
is not physically justified and it introduces an addi- 
tional spurious dependence on the parameter q into 
the problem. 

Let us briefiy summarize the consistent system of equa- 
tions. The energy density wb of a magnetic field of inten- 
sity B per volume unity is 



where /n is a geometrical factor of the order of unity. If 
the magnetic field produces some instability with a ver- 
tical component, one must have ub > un- Otherwise, 
the restoring force of gravity which acts at the dynamical 
timescale would immediately counteract the magnetic in- 
stability. From this inequality, one obtains P < 



2/n 



Ar2 • 



If /n = 2 ' have the condition for the vertical amplitude 
of the instability f Spruit .2002. Eq. 6), 



I < r 



UJA 

N 



(5) 



where r is the radius. This means that there is a maximum 
size of the vertical length / of a magnetic instability. In 
order not to be quickly damped by magnetic diffusivity, 
the vertical length scale of the instability must satisfy 



CTB 



77 n 



(6) 



where fl is the angular velocity and ctb is the characteristic 
frequency of the magnetic field. In a rotating star, this 
frequency is ctb = (w^/fi) due to the Coriolis force (Spruit 
nUU^ see also Pitts & TaylerOWSI). The combination of 
the limits given by Eqs. (0) and © gives for the case of 
marginal stability. 



1 2 2 



with 



B 



UJA 



(4 np) 2 



(1) 



where ua is the Alfven frequency in a spherical geometry, 
r is the radius and p the density. In stable radiative lay- 
ers, there is in principle no particular motion. However, 
if due to the magnetic field or rotation, some unstable 
displacements of vertical amplitude 1/2 occur around an 
average stable position, the restoring buoyancy force pro- 
duces vertical oscillations around the equilibrium position 



n 



(7) 



The equality of the amplification time of Tayler insta- 
bility Ta = N/{u>a^q) with the inverse of the charac- 
teristic frequency ctb of the magnetic field leads to the 
equation (Spruit EOOSl 



n 



N 



with 



d\nr 



(8) 
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Fig. 1. Internal distribution of the angular velocity fl{r) 
as a function of the lagrangian mass in solar units in a 
15 Mq model, without magnetic fields, at various stages 
of the model evolution indicated by the central H-content 
Xc during the MS-phase. The initial velocity Ujni = 300 
km s~^. 



With account of the Brunt- Vaisala (Eq. one has 



^^T r,/K 



(9) 



By eliminating the expression of N'^ between Eqs. ^ and 
we obtain an expression for the magnetic diffusivity. 



(10) 



Eqs. 10 and form a coupled system relating the two 
unknown quantities r; and loa- Instead, one may also con- 
sider for example the system formed by Eqs. lO and H10() . 
Formally, if one accounts for the complete expressions of 
the thermal gradient V, the system of equations would be 
of degree 10 in the unknown quantity x — (^) (paper 
II). The fact that the ratio rj/K is very small allows us to 
bring these coupled equations to a system of degree 4, 



2n^q^ = . (11) 



The solution of this equation, which is easily obtained nu- 
merically, provides the Alfvcn frequency and by Eq. (|10|l 
the thermal diffusivity. As shown in paper II, the various 
peculiar cases studied by Spruit r2002|l are all contained 
in the more general solution given here. 
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Fig. 2. Same as in Fig. ^ but with magnetic fields. One 
notices the almost constant values of the angular velocity 
in the models. 



3. Coefficients for the transport of angular 
momentum and chemical elements 

The condition that the ratio loa/VI given by (Eq. ISJ is 
equal to or larger than the minimum value defined by Eq. 
(|7|) leads to a condition on the minimum shear for the 
dynamo to work (cf. Spruit | 



1/4 



(12) 



There, N is given by Eq. (0) above and rj by Eq. l(Tn|l . If 
this condition is not fulfilled, there is no stationary sit- 
uation and we consider that there is no magnetic field 
present. In practice, we will see that this situation occurs 
in the outer stellar envelope. There is a second condition. 
We need to check that ft > uja- If this last condition is 
not realized, the present system of equations does not ap- 
ply and we should consider the case of very slow rotation. 
This case is treated in the Appendix. In the present mod- 
els, the value of considered is large enough so that we 
did not have to apply the solution for the case of very slow 
rotation. 

The azimuthal component of the magnetic field is 
much stronger that the radial one in the Tayler-Spruit 
dynamo. We have for these components (Spruit P002) 



(47rp)2 r UJA and B^- = (IJr) , 



(13) 



where uja is the solution of the general equation (|lll) and 



Ir is given by 



N 



which is obtained by assuming 



the marginal stability in Eq. ((SJ. 
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Turning towards the transport of angular momentum 
by magnetic field, we first write the azimuthal stress by 
volume unity due to the magnetic field 




Then, the viscosity v for the vertical transport of angular 
momentum can be expressed in terms of S (Spruit |2002), 

This is the general expression of v with ljx given by the so- 
lution of Eq. ((nil and with N by Eq. Q. We have the full 
set of expressions necessary to obtain the Alfven frequency 
WA and the magnetic diffusivity 77. The parameter 77 also 
expresses the vertical transport of the chemical elements, 
while the viscosity v determines the vertical transport of 
the angular momentum by the magnetic field. 

In paper II, we also checked that the rate of magnetic 

energy production — J^'W ~ iP'^^'^'W P*'^ '^^ 
time and volume is equal to the rate — ^pvVl^q^ of the 
dissipation of rotational energy by the magnetic viscosity 
V as given above. 

4. Numerical applications 

4.1. The models 

We consider here models of 15 Mq, with a standard com- 
position of X = 0.705 and Z = 0.02. The physics of the 
models, opacities, nuclear reactions, mass loss rates, struc- 
tural rotational effects, shear mixing, meridional circula- 
tion, etc. are the same as in recent models of the Geneva 
group (Meynet & Maeder IW)^ . We calculate 3 sets of 
models: the first is without rotation, the second with ro- 
tation and the third one with both rotation and magnetic 
field. The initial rotation velocity is 300 km s~^, which 
leads to an average velocity on the MS phase of about 240 
km s^^ when no magnetic field is present (see Fig. |3J)- In 
paper II, we did not include the effects of the meridional 
currents, in order to independently study the effects of the 
magnetic field. However, we showed that meridional cir- 
culation and the magnetic field may significantly interact, 
with possible consequences for the transport mechanisms. 
Thus, to examine this interplay, we also account here for 
meridional circulation currents. For the moment, we do 
not account for magnetic coupling by stellar winds, since 
there is no external convective zone in massive stars on 
the Main Sequence. 

The choice of the time step is imposed by the most 
rapid process taking place. We checked that changing the 
time steps did not change the results. Here the fastest 
process is the transport of the angular momentum by the 
magnetic diffusivity. The diffusion coefficients v for the 
angular momentum by the magnetic field reach in some 
cases values of 10^^ cm^ g~^, while the average value is 1-2 
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Fig. 3. Evolution of the rotation velocities at the surface 
of 15 M0 models during the MS phase with and without 
magnetic fields H. The initial velocity Uini — 300 km s~^ 
in both cases. We see the much higher surface rotation 
when magnetic field is included. 

orders of magnitude lower. The cases of strong coefficients 
lead to diffusion timcscales r « (A_R)^/i/ « a few 10^ yr, 
where D is the diffusion coefficient. For an appropriate 
treatment, we need to adopt very small time steps. In 
practice, we take time steps of the order of 20 years. This 
implies about 6 • 10^ models to cover in a detailed way the 
exact interplay of the effects of the magnetic field and of 
the meridional circulation during the MS phase ! 

Clearly, the model properties do not change signifi- 
cantly over such a time scale as shown for example by 
Fig. 121 Thus, in future, faster processes of calculations 
may be devised. However, by imposing very short time 
steps on the two processes, we set the study in the linear 
regime for both instabilities and we may thus proceed to 
an addition of their own different effects. 

4.2. Evolution of the internal and superficial rotation 

In a rotating star, the internal profile of angular velocity 
n{Mj-) changes with time during evolution, due to vari- 
ous effects (Meynet & Maeder JOOOJ: central contraction 
and envelope expansion, transport of angular momentum 
by circulation and mass loss at the stellar surface. One 
assumes here that the mass lost by stellar winds just em- 
barks its own angular momentum. The case of anisotropic 
stellar winds has been studied by Maeder H2002|l . who has 
shown that such effects are important only for very mas- 
sive stars {M > 60 M0) in fast rotation. Here, the loss of 
angular momentum at the stellar surface only has a lim- 
ited importance for the evolution of rotation. Fig.Qshows 
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Fig. 4. Distribution of the magnetic field in the model of Fig. 5. The diffusion coefficients in the model with central 
initial 15 Mq when the central hydrogen content is Xc = H-content Xc = 0.60. In the outer envelope, where the 
0.60. The age of th model is 4.2798 • 10^ yr and the actual magnetic field is absent, the transports by shear and by 



mass is 14.97 Mp 



0- 



combination of the meridional circulation and horizontal 
turbulence expressed by Doff are present. 



the evolution of the internal profile of rotation. Differential 
rotation builds up during the MS phase, reaching about 
a factor of two between at the surface and in the con- 
vective core, near the end of the MS when the central 
H-content is Xc = 0.10. Then, fast central contraction 
accelerates the core rotation. 

The situation is quite different when magnetic fields 
as described above are accounted for. As shown by Fig. 
12 the angular velocity fl is almost constant throughout 
the stellar interior. It is not exactly constant, otherwise q 
would be zero and the magnetic field would not be sus- 
tained anymore. Only at the very end of the MS phase 
does the fast central contraction bring about a small sig- 
nificant difference. 

Fig.|31shows the evolution of the surface rotation veloc- 
ities as a function of age in the models with and without 
magnetic field. We notice the higher velocity during the 
whole MS phase of the model with magnetic field, which 
has V « 300 km s~^, compared to the model without a 
field, where v « 250 km s^^. There are 3 effects interven- 
ing: 

1. The main part of the difference already occurs on the 
ZAMS. With rotation only, the model adjusts its rota- 
tion very quickly to an equilibrium profile determined 
by meridional circulation (Meynet & Maeder l!jOOO|l and 
this establishes an fi-gradient which reduces the sur- 
face velocity. On the contrary, with a magnetic field, 
the model can keep a constant Q and there is no initial 
decrease. 



2. The other two smaller effects occur during MS evolu- 
tion. The magnetic coupling forces the surface to co- 
rotate with the convective core, as contraction makes 
it spin-up. With rotation only, the coupling (due to 
meridional circulation) is much weaker. 

3. In the magnetic model, the surface enrichment in he- 
lium is slightly higher, this keeps the radius smaller. 
Thus, as the radius expansion is smaller, the decrease 
of surface rotation is also smaller. 

The question of the differential rotation in the hori- 
zontal direction has also to be examined. In current ro- 
tating models without a magnetic field, it is assumed that 
the horizontal turbulence is strong enough to suppress the 
horizontal differential rotation, so that the hypothesis of 
shellular rotation with fl = fl{r) applies fZahn [T?)92|l . In 
the presence of a magnetic field, the horizontal turbulence 
is likely reduced. The numerical values for the field 
found below is of the order of a few 10'' G (cf. Fig.^J. The 
horizontal coupling ensured by the magnetic field (cf. Fig. 
Q is expressed by a coefficient Z^Bh (Maeder & Meynet 

Db,^ - r2 {uji/n) = — ^ for Q ^ uja . (16) 

The value of Db,, is of the order of 10^^ cm^/s in the 
present 15 Mq model. Thus, the magnetic field is large 
enough to ensure the horizontal coupling, so that the as- 
sumption of shellular rotation is still valid. 
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Fig. 6. Internal distribution of the hydrogen mass frac- 
tion X as a function of the Lagrangian mass in a 15 M0 
model with rotation at various stages of the model evo- 
lution from top to bottom during the MS-phase. Some 
convective zones appear in the H-shell burning stage at 
the end of the MS phase. 



Fig. 7. Internal distribution of the hydrogen mass fraction 
X as a function of the Lagrangian mass in the model with 
magnetic field calculated according to expressions of the 
present paper at various stages of the model evolution 
from top to bottom during the MS-phase. 



4.3. Diffusion coefficients for the transports of the 
angular momentum and chemical elements 

We need to determine in which regions of the star the 
magnetic field is present, its intensity and the run of the 
various diffusion coefficients. As an example, we examine 
the 15 Mq model at the beginning of the evolution, when 
the central H-content is Xc — 0.60. In this model as in 
all MS models, q > qmin in a region which starts just 
above the core and extends over most of the envelope. In 
the model with Xc — 0.60, the magnetic field is present 
from A/j. = 6.7 Mq to 13.4 Mq, as illustrated in Fig. H 
The average field is about 2 • 10"* G and the corresponding 
value of ojA/n is about 5- 10~"'. Above 13.4 Mq, condition 
H12I) is no longer satisfied. The region where magnetic field 
is present slightly increases during MS evolution as the 
convective core recedes. Near the end of the MS phase, 
when Xc = 0.05, the magnetic field is present from 4.2 
Mq up to 13.8 Mq. 

The runs of the various diffusion coefficients are il- 
lustrated in Fig. [S| The largest diffusion coefficient is i' 
which acts for the vertical transport of angular momen- 
tum, the large value of v imposes the nearly constant Q in 
the interior. The value of is about 6 orders of magnitude 
larger than the diffusion coefficient 77 for the transport of 
the chemical elements, so that the surface enrichments 
in products of the CNO cycle due to this magnetic dif- 
fusion coefficient only are rather limited. The coefficient 
Doff which applies to the transport of chemical elements 



by meridional circulation is much larger than -Dshoar, while 
in rotating stars without magnetic field, it is generally the 
opposite. This is due to the fact that the velocity of merid- 
ional circulation is much larger when O is almost constant 
throughout the stellar interior (about 10~^ cm s~^), while 
-Dshoar is much smaller. Indeed Dghear is about 4 orders of 
magnitude smaller than in rotating stars without a mag- 
netic field, as a consequence of the near solid body rotation 
in magnetic models. 

Fig. |S1 also shows that the value rj/K is always very 
small, which justifies the simplifications made in deriv- 
ing Eq. I|ll|l . Apart from the slight extension of the zone 
covered by the magnetic field as mentioned above, the 
orders of magnitude and relative ratios of the diffusion 
coefficients remain about the same during the whole MS 
evolution. 

An important result of Fig. is that the transport by 
meridional circulation (expressed by Dcs) is in general 2-3 
orders of magnitude larger than the transport of the ele- 
ments by the magnetic field (expressed by r/) . This clearly 
shows that the direct effect of the magnetic instability is of 
little importance with respect to thermal instability, which 
drives meridional circulation, for the transport of chemical 
elements. In addition, meridional circulation is important 
for the transport of the elements from the external edge of 
the magnetic zone to the stellar surface. However, for the 
coupling of Q, the magnetic field is much more efficient 
than meridional circulation. One can define a velocity Ub 
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Fig. 8. Top lines near 15 M©: evolution of the stellar mass. 
Lower lines: evolution of the mass of the convective core. 
The dotted line applies to the model without rotation, 
the short-broken line to the model with rotation {vini— 
300 km s^^) but without magnetic fields, the continuous 
line to the model with rotation (wini= 300 km s^^) and 
magnetic fields. 

for the vertical transport of angular momentum by the 
magnetic field (cf. paper I) 

u d\nil , , 

Ub = 5- —— . (17) 
r omr 

This quantity may be compared to the vertical component 
of the velocity U{r) of meridional circulation. In general 
Ub is much larger than U{r), which confirms the dominant 
role of the magnetic field for the transport of angular mo- 
mentum. 

4.4. Internal distribution of hydrogen 

Figs. El and 13 show the internal distribution of hydrogen in 
the cases without and with magnetic field. Firstly, rotation 
with a magnetic field leads to the formation of a larger 
core at the end of the MS phase compared to models with 
rotation only (models with rotation themselves have larger 
cores compared to models without rotation, cf. also Fig. 

In models with rotation only, there is a significant ero- 
sion of the /i-gradient at the immediate edge of the core, 
developing during the second half of the MS phase. This 
erosion directly at the edge of the core is absent in models 
with a magnetic field. The reason is Eq. (|10(l . which gives, 
if the /i-term dominates in the Brunt- Vaisala frequency , 
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Fig. 9. Evolutionary tracks in the HR diagram for 15 
M0 stellar models. The dotted line applies to the model 
without rotation, the short-broken line to the model with 
rotation (wini= 300 km s^^) but without magnetic fields, 
the continuous line to the model with rotation (uini= 300 



km 



and magnetic fields. 



(18) 



As noted in paper II, the mixing of chemical elements de- 
creases strongly if the /i-gradient grows and this effect lim- 
its the chemical mixing of elements by the Tayler-Spruit 
dynamo in the regions just above the convective core. 

Finally, in models with a magnetic field the mixing in 
the envelope is much greater than in the case without a 
field. This effect appears during the entire MS phase and 
the final helium surface content reaches 1^ = 0.31. This 
effect is not due to the shear, which is essentially absent 
in magnetic models. It is also not due to the magnetic dif- 
fusivity, because it is rather small and in the very external 
regions there is no magnetic field. However, the transport 
of chemical elements by the circulation (an effect described 
by Dcb) is larger and is generally the main effect in the 
envelope. It produces a slight transport of the elements, 
which enriches the stellar surface in elements of the CNO 
cycle (cf. Fig. llOfl . The models of paper II with a magnetic 
field, but without meridional circulation predict no N ex- 
cesses and this was a difficulty. However, magnetic models, 
when the meridional circulation is included, lead to a sig- 
nificant N-enrichment, which compares with observations 
(cf. Sect. 123. 



5. Evolution of observable parameters 

The magnetic field created by the dynamo does not reach 
the stellar surface. However, there are consequences con- 
cerning the tracks in the HR diagram and especially the 
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Fig. 10. Time evolution of the surface helium content Yg in mass fraction, of the N/0, N/H and N/C in mass fraction 
for various models: The dotted line applies to the model without rotation, the short-broken line to the model with 
rotation {vini— 300 km s^^) but without magnetic fields, the continuous line to the model with rotation (uini= 300 
km s^^) and magnetic fields. 



enrichment of surface abundances, which may lead to ob- 
servable consequences. 

5.1. Evolution of the core, lifetimes and tracks in the 
HR diagram 

The evolution of the tracks in the HR diagram is deter- 
mined by the evolution of the core mass fraction. Fig. |S1 
shows the evolution of the total stellar mass in the mag- 
netic case. When — 0.30, the actual mass is 14.875 Mq. 
Mass loss increases near the end of the MS phase and the 
mass is 14.389 Mq when = 0.000. As is well known, 
rotational effects slightly increase the core mass fraction 
near the end of the MS phase and produces longer MS 
lifetimes. The account of the magnetic field further en- 
larges the core mass fraction as well as the MS lifetimes 
as illustrated in Fig. |H1 

The effects on the tracks in the HR diagram are shown 
in Fig. El In addition to the structural changes, rotation 



also produces a distortion of the stellar surface and an in- 
crease of the average radius (estimated at P2(cosi9) — 0, 
where is the colatitude). This makes a red wards shift 
of the tracks of rotating models in the HR diagram, espe- 
cially visible near the ZAMS, since no other effects inter- 
vene there. For the model with rotation only, we see the 
upwards shift due to the slightly larger core and a slightly 
bluer track resulting from a tiny surface He-enrichment. 
The model with rotation and magnetic field reaches a 
higher luminosity at the end of the MS due to the larger 
core. It is also bluer due to the larger surface enrichment in 
helium, which lowers the opacity and decreases the stellar 
radius. It is doubtful that the above effects are sufficient 
to infer the presence or absence of magnetic fields. 
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5.2. Evolution of the abundances of helium and CNO 
elements at the stellar surface 

The evolution of the surface abundances in helium and 
CNO elements results from the internal profiles shown in 
Figs.Eland[71 Fig.llOlshows for the 3 models considered the 
evolution of the helium abundance in mass fraction at 
the surface, the N/H, N/0 and N/C ratios of numbers of 
atoms. We see the same trends in the four panels. There is 
no enrichment in absence of rotation. With rotation only, 
there are moderate enrichments, by a factor of 2 for N/H, 
2.5 for N/C and 1.8 in N/0. With both rotation and 
magnetic field, the surface enrichments are much larger, 
particularly for the helium abundance which reaches Yg 
0.31 at the end of the MS phase. The increases in N/H 
reach a factor of 5, 11 for N/C and 6 for N/0. 

Observationally, there are many estimates of N en- 
richments for OB stars. A recent review of the subject 
has been made by Herrero & Dufton H2()()4(l . who show 
clear evidence of rotationally-induced mixing in OB stars, 
also with the result that some fast rotating stars do 
not show N~enrichments. Four 09 stars were studied by 
Villamariz et al. H2002(l : three low rotators have an excess 
of [N/O] < 0.2 dex, one fast rotating 09III star with 
vsini = 430 km has an excess of [N/O] = 0.7 dex. 
The two most massive stars with fast rotation in the asso- 
ciation Cep 2 have an excess of [N/O] < 0.3 dex ( Daflon 
et al. 12001(1 . Two stars in Sher 25 have excesses in [N/O] 
by a factor of 3 to 4 (Smartt et al. 12002(1 . Observations of 
the N/H ratios by Venn and Przybilla 1(2003(1 for galac- 
tic A-F supergiants show an average excess of a factor of 
3, with extreme values up to a factor of 8. The orders of 
magnitude of the predicted and observed enrichments are 
similar. However, the situation is still uncertain due to the 
relative lack of accurate observational data for MS stars. 



6. Conclusions 

The main result is that a magnetic field imposes nearly 
solid body rotation and this favours higher rotational ve- 
locities during MS evolution compared to cases where the 
magnetic field is not accounted for. Due to the nearly con- 
stant in the stellar interior, the transport of chemical 
elements by shear mixing is negligible. The transport of 
elements by Tayler-instability is also very limited. An in- 
teresting feature of the model is that meridional circula- 
tion is strongly enhanced by the flat f2-curve and this is 
the main effect influencing the transport of the chemical 
elements in the present models. 

There remain however some doubts as to whether the 
Tayler-Spruit dynamo is really active in stellar interiors. 
Up to now, magnetohydrodynamic models expressing the 
growth and evolution of the magnetic field in rotating 
stars have not yet confirmed the existence and efficiency 
of this particular instability in stellar interiors (Mathis 
12005(1 ■ Even in the observationally and theoretically much 
better studied case of the Sun, the exact location, origin 



and evolution of the solar dynamo are still not fully un- 
derstood. 

The larger size of the predicted surface enrichments is 
well (or even better) supported by observations. However, 
this remains uncertain in view of the small number of ac- 
curate observations. In this respect, it might be crucial 
to also have such comparisons for stars in the Magellanic 
Clouds where the observed enrichments are much larger 
than in the Galaxy. The answer may come from astero- 
seismology. There are no p-modes expected, however g- 
modes may be present and yield some information on the 
internal r2-distribution. At present, this seems to be the 
most compelling possible test. 
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Appendix A: Equations for the case of low 
rotation 

The dynamo properties were established (Spruit EES and 
Sect. 2) for the case of fast rotation, with the condition 
that the rotation rate is larger than the Alfven frequency 
LUA, i.e. Q 3> UJA- In the fast rotating case, the growth rate 
of the magnetic field is reduced by a factor lua/^ as firstly 
suggested by Pitts & Tayler 1(1986(1 . The above equations 
of the dynamo have been derived under this condition. 
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However, if we want to study very slowly rotating stars, 
we need also to consider the case where 



(A.l) 



This may be interesting also in central regions of radia- 
tive stars. There, as a result of the small value of r, the 
Alfven frequency uja becomes large so that the regime cor- 
responding to Eq. I)A.1(I may apply instead of the usual one 
(cf. Sect. EI). 

The energy of the Tayler instability fTavler [l973|l must 
be large enough to overcome the restoring force of buoy- 
ancy and this implies that the Alfven frequency must be 
larger than some limit depending on N. The vertical ex- 
tent Zj. of the magnetic instability is also given by Eq. (O 
above. In order that the perturbation is not too quickly 
damped by the diffusion of the magnetic field, one must 
have > — . At low rotation, in the absence of significant 
Coriolis force, the frequency ctb = ^A- The combination 
of these two limits in the low rotation case leads to the 



condition i 



> iV- 



2 V 



and for the marginal situation cor- 



responding to equality, we have 



r2 n 



(A.2) 



This equation relates the magnetic diffusivity rj and the 
Alfven frequency uja instead of Eq. |(2J . 

The field amplitude may be fixed by the equality of 
the amplification time Ta of the Tayler instability and of 



the timescale cTg of the field (Spruit 



One has 



(A.3) 



Instead of the expression in ((HJ, we get for the low rotation 
case with ctb = 



N 



(A.4) 



UJA has disappeared from the equation. How should we 
interpret this ? This can be done with the help of Eq. 10) , 
which gives 



rj/K + 2 n 



V 2 



(A.5) 



The two equations (|A.2|) and (|A.5|) form our basic system 
of equations with the two unknown quantities ry and uja- In 
order to have a diffusivity rj positively defined, Eq. (jA.Sp 
imposes that the differential rotation parameter q must be 
larger than some minimum value Qmim 



(A.6) 



If this is the case, then according to Eq. (|A.2|) the same is 
true for uja and thus a magnetic field is effectively created 
by the slow dynamo process. 



Let us now estimate the transport coefficients. 
Eq. (|A.5p immediately leads to the following expression 
for the magnetic diffusivity 77, 

(A.7) 



V 



In all cases, the ratio rj/K is very small, typically of the or- 
der of 10~^ (cf. Fig.O. Thus, we obtain the corresponding 
expression for the diffusivity. 



V 



2 K 



72 



(A. 



This equation provides the magnetic diffusivity rj, if \q\ is 
larger than |gmin|- The coeflacient 77 determines the trans- 
port of the chemical elements by the magnetic instability. 

Now, with Eq. l|A.2p . we obtain the Alfven frequency 
UJA at each location r in the star, 



ujaY _2K 



r2 n 



iVZ^2 



(A.9) 



This confirms that for a star rotating with angular velocity 
51 and having a certain /i-gradient, a magnetic field is 
created only if the actual value of the differential rotation 
parameter \q\ is larger than Igminl- 

The azimuthal stress S due to the magnetic field gen- 
erated by Tayler instability becomes 



S — - — B^B^p — 

4 TT 



1 

4 TT 



p r 



N 



(A.IO) 



The viscosity v is expressed in terms of 5* and we get 

S _ nr^ /ujaY 
pqn ^ \N 



V = 



(A.ll) 



With Eq. (jA.9|) for the Alfven frequency, we get finally 



2 K 



(A.12) 



which determines the transport of angular momentum. 
The expression of v is the same as the expression of 77 
fEa. lA.8|) . This means that in the case of low rotation, the 
Tayler magnetic instability transports in a similar way the 
chemical elements and the angular momentum, while at 
high rotation the magnetic coupling is stronger by a fac- 
tor than the transport of the elements, according 

to the general expressions by Maeder & Meynet H2l)U4f) . 
This factor comes from the expression of ctb, and as in 
the present case the ratio (^^) is absent in ctb, wc get 
equality of v and rj. 

Condition (|A.4p must also be satisfied in order to have 
angular momentum transport. The equations (jA.7|) or 
(IA.8|) and (|A.12I) provide the transport coeffcients at each 
stellar layer as a function of the local quantities, such as 
57, g, N'^, N'^, etc... available in the stellar models. The 
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Alfven frequency and the magnetic field intensity are ob- 
tained by Eqs. (|A.9|I and (jlj. 

Tlie rate of magnetic energy production per unit 
of time and volume must be equal to the rate of the 
dissipation of rotational energy by the magnetic viscosity 
V as given above. This check of consistency was verified 
for the case of high rotation. One assumes here for sim- 
plification that all the energy dissipated is converted to 
magnetic energy. One has 

= ^pi^n^q^ 

which with Eq. IIA.12p gives the dissipation rate, 

w.^pn\'^{q'n'-N'^). (A.14) 

The magnetic energy per unit volume is , it is pro- 
duced in a characteristic time given by = ^a^- Thus, 
one has 

Wb = ^iOA = Ipr'ui = P1^^'§2 (.1' - ' (A- 15) 

where we have used the Eq. (|A.9|I for the Alfven frequency. 
Thus, the expressions for Wi, and Wb are the same. We 
see that no magnetic energy is produced if condition (jA.6|l 
is not realized. This shows the consistency of the field 
expression for B^p and of the transport coefficient v. 



(A.13) 



